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Abstract 

We derive strong contraints on the Yukawa couplings and the vacuum expectation value 
in the singlet majoron model, taking into account the possibility of a small gravitationally 
induced mass for the majoron. If the present baryon asymmetry was created earlier than 
at the electroweak phase transition, then to preserve it from the combined effects of lepton 
number violating interactions and electroweak sphalerons, we find that stringent constraints 
on the Yukawa couplings, h < 10~ 7 , and the majoron scale, v s < vew, must be satisfied. We 
also carefully rederive baryogenesis bounds on neutrino masses, finding that in general they 



apply not to the masses themselves, but only to related parameters, and they are numerically 
somewhat less stringent than has previously been claimed. 



1 Introduction 



There is an ever increasing interest in the possibility for neutrino masses. Aside from the 
intrinsic interest in new and perhaps tangible physics, neutrino masses have been shown 
to be quite useful. A neutrino with a small mass, O(10 eV), could make a significant 
contribution to the overall energy density of the Universe. Other potential benefits from 
neutrino masses include MSW neutrino mixing relating to the solar neutrino prolem [|IJ , and 
the possibility of exploiting lepton number violating interactions in conjunction with baryon 
and lepton number violation in the electroweak sector to generate a baryon asymmetry 
of the Universe |3j . 

The very small ratio of expected neutrino masses, relative to other particle masses, made 
the idea of the seesaw mechanism quite attractive Numerous models incorporate a 
neutrino mass seesaw, including unification via SO(10) ||, or simply adding a single complex 
singlet to the standard model which breaks a global XJ(1) b-l symmetry, ie., the singlet 
majoron model ||. In both of these models lepton number is spontaneouly broken. In a 
GUT, {B — L)-violating effects will generally be small because the breaking scale is generically 
Mqut- In the majoron model, there are a priori few restrictions on the breaking scale, which 
might be close to the electroweak scale. 

However the possibility of lepton number violation at or above the electroweak scale is 
potentially dangerous, as was first discussed by Fukugita and Yanagida in the context of 
a simple model with explicit lepton number violation. The danger occurs when L-violating 
interactions are in equilibrium in the early Universe at the same time that (B + L)-violating 
nonperturbative electroweak interactions, due to sphalerons, are in equilibrium. In this 
case the baryon asymmetry of the Universe is destroyed independently of the initial B — L 
asymmetry. This realization has led to a series of contraints on models with additional 



sources of baryon or lepton number violation 0-[l3]. These bounds are of course evaded if the 
baryon asymmetry is ultimately produced during or after the electroweak phase transition. 

The baryon wash-out bounds also directly apply to the majoron model if the lepton 
symmetry breaking occurs above the electroweak scale. They take on an especially interesting 
twist however, when combined with the "expectation" of gravitational corrections to the 
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singlet sector. Recalling the well-known adage that Nature abhors global symmetries, it is 
quite reasonable to expect nonrenormalizable contributions to the low energy Lagrangian due 
to gravitational corrections that violate nongauged symmetries jnj. Although such terms 
explicitly break global symmetries, their effects are strongly felt only above some cutoff, in 
this case the Planck scale. Such effects were previously considered in the triplet majoron 
model |L5] and for axions |16| . Nonrenormalizable gravitationally- induced interactions which 
explicitly break the global B — L symmetry were recently studied in the context of the singlet 



majoron model |T7|, [18| . In [|Tj] it was shown that these effects generate a mass 0(1-2) keV 
or larger for the majoron and that the cosmological consequences of such a massive majoron 
imply a stringent upper bound on the vacuum expectation value of the singlet field in the 
broken phase. 

In this paper we derive some of the baryon wash-out bounds in detail. In particular 
we extend the bound coming from decays and inverse decays Jll| of heavy neutrinos to the 
case of spontaneously broken symmetry and find an explicit expression for this bound in 
terms of the light neutrino seesaw masses, the neutrino mixing matrix and the initial lepton 
asymmetries. We will then show how the bound on the singlet vacuum expectation value 
may be further strengthened when combined with the (B + L) -violating nonperturbative 
electroweak effects. 



2 The Model 

The singlet majoron model extends the Standard Model by merely adding a single complex 
scalar field, cr, and a two-component fermion field ur, both SU(2)xU(l) singlets. The new 
interaction Lagrangian terms may be written as 

hHLv R + h 2 av R VR + h.c. + V(a, H); (1) 

V(a, H) = X s (\a\ 2 - v 2 s /2) 2 + \ H (\H\ 2 - v 2 /2) 2 + X SH (\a\ 2 - v 2 s /2){\H\ 2 ~ v 2 /2), (2) 

where we allow for self-interactions of the scalar singlet as well as interactions with the 
Standard Model Higgs doublet, H. L is the usual left-handed lepton doublet, with generation 
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indices suppressed. When a aquires a vacuum expectation value we can write 

v = -j={v* + P + ix)- (3) 

Upon symmetry breaking, vr acquires a Majorana mass, M = \[2h2V s and a Dirac mass 
with vl, m = hv/y2 where vj a/2 = 174 GeV is the VEV of H. For M m the low energy 
mass eigenstates are determined by the seesaw relations iV = vr + i>r c + (m / 'M){i>l, c + vl) 
and v = ul 4~ v^ c — {m/M)(vR c + ur) with m^v — M and ~ m 2 /M. The field p has mass 
m p ~ f s for self-couplings of order 1, as is often assumed, and the massless majoron is \. 

Because this model introduces only new singlet fields, it is only weakly constrained. If any 
of the light neutrinos are heavier than ~ 100 eV, they must decay fast enough so that neither 
they nor their decay products overdose the Universe fT9| . Furthermore, they must decay fast 



enough to allow for sufficiently long late-time matter domination, so that primordial density 



fluctuations can grow to form galaxies [20|. The observed neutrino pulse from SN 1987A 



leads to a constraint on the VEV, v s for neutrinos in the mass range from 100 eV to ~ 30 



MeV [21 1 . A massless majoron would also conflict with bounds on light particle degrees of 



freedom [22| unless it decoupled early enough f23|. For the very small neutrino masses we 



will argue for below, the majoron was always sufficiently decoupled. 



3 The Gravitationally Induced Bound 



The gravitationally-induced interactions for the singlet majoron model, considered in ||17|| . 
were assumed to have the generic form 

with the consequence that for v s > v, a majoron mass 

3/2 



V20X (^A keV (5) 



is generated. Similarly, although we will not be interested in this case, when v s < v, a 
majoron mass of at least O(keV) will also be produced. We anticipate that A ~ 1 (there is 
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no reason to assume otherwise), but will nevertheless keep all A-dependence explicit in our 
results. 

An immediate consequence of majorons with a mass in the keV range is that they would 
overdose the universe, if they were stable. To see this, one must estimate their relic abun- 
dance, Y = n x /n 1 , given by the ratio of particle species in equilibrium today to the number 
at the time they decoupled. Above the phase transition where a gets its VEV, all the parti- 
cles are massless, so that majorons are in equilibrium at least until the critical temperature 
T c , which is of order v s . If this occurs above the weak scale, but below the threshold for new 
SUSY or GUT particles, then the spectrum of particles in equilibrium will just be all those 
of the Standard Model, resulting in the dilution of the majoron number density by a factor 
Y ~ (1/27) p2[ . This dilution is sufficient to satisfy the constraints from big bang nucleo- 
synthesis |24[] since the the majoron at the time of nucleosynthesis would only contribute 



the equivalent of ~ 0.03 neutrinos [p2 , well below the bound 0.3. However, the present 
density of majorons becomes 

Px = 2 m x Yn y ( 6 ) 



m x < 340 eV (-U (7) 



which leads to the constraint 



\27Yj \0.25 / 

and which for A £ 0.01 (v/v s ) 3 is incompatible with the estimate in (||) [|17"| . If the majorons 
are kept in equilibrium to lower temperatures, for example through the processes N ^ x v 
or vv <-> xXi Y wn l b e larger and the conflict will only be exacerbated. Note that the bound 
cannot be evaded by assuming the majoron decoupled at a higher temperature (resulting in 
a smaller value of Y), since the majoron mass grows faster as a function of v s ~ T c than does 
Y~ l . Thus majorons must decay. But not only must they avoid overdosing the universe, 



they must decay early enough to allow for the growth of density fluctuations [20 



The recent measurements by COBE [25| of the quadrupole moment in the spectrum of 



primordial density fluctuations, which is consistent with a flat spectrum and an amplitude 
y ~ 5 x 10~ 6 , indicates the need for a substantial period of matter domination to allow for 
the growth of these perturbations in order to form galaxies. Because of the majoron mass 
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scale, m x Y £0.1 keV, galaxy size perturbations will go nonlinear only if the last epoch of 
matter domination is sufficiently long [fit], With such a small value of ^ it is necessary 
to assume that majorons never dominate the overall density. This implies that [2D| 

^ £ 9.0 (8) 

where To is the majoron decay temperature. Relating To to the majoron lifetime gives the 
limit 



m v Y \ f t d \ 1 / 2 4 



1 keV/ \1 sec 



< 1 x 10 4 . (9) 



Then using 



TD =rD= 16^ I" ) (10) 
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one can obtain the bound [17] on v s , 



v v < 20 TeV — A" 1 / 7 (11) 

1 25 eV/ V27F/ V ' 



which is interesting in and of itself (note the weak dependence on A), as long as m v < 25 eV. 
In fact the latter condition must be satisfied, because m v here refers to the heaviest neutrino 
which can result from the decay x ~^ vv - Such a neutrino, if heavier than the cosmological 
bound of 25 eV, would itself have to decay into three lighter neutrinos, v — > 3z/, with a 
majoron in the intermediate state. However the vv'x coupling is highly suppressed [p7 [; at 
low energies the Yukawa interactions of the singlet fields have the form 

fr(*xs-&x, + *,x, + o(£),x,), ( 12) 

where X = p + ij^X- Because the rate is suppressed by the tiny factor (m/M) 12 = (m u /M) 6 , 
this decay mode is far too slow to vitiate the 25 eV bound on m v . 



4 Sphaleron Induced Bounds 

We now turn to the baryon wash-out constraints. Any model with lepton number violation 
faces the danger of baryon number wash-out as a result of the combination of its own 
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lepton number violating interactions and the (B + L) -violating sphaleron interactions in 
the Standard Model. This was first noticed by Fukugita and Yanagida p], who showed that 



the presence of a lepton (or baryon) number violating operator such as j^vvHH will erase 
any primordial baryon asymmetry if these interactions are in equilibrium simultaneously 
with sphaleron interactions. The condition that these interactions be out of equilibrium 
translates into a bound on the neutrino mass [0, || |9[ 

m „< 50k ev(l^X) 1/2 (13) 

where Tbl is the temperature at which a B — L asymmetry was produced or ~ 10 12 GeV 
(the maximum temperature for which sphaleron interactions are in equilibrium), whichever 
is smaller. This type of bound has been extended to general baryon and lepton number 



violating operators in |10| where some general constraints on the majoron model were also 
considered. 

Another, complementary bound on neutrino masses can be obtained considering the 
lepton number violating decays and inverse decays iV «-> HL, that inevitably accompany any 
nonrenormalizable operator such as indicated above. The requirement that these processes 
are out of equilibrium has been shown to imply |T| 



m v < 1(T 3 eV (14) 



if the B — L asymmetry is generated above T ~ M. Remember that both fll3|) and ([14]) can 
be evaded if the baryon asymmetry is ultimately produced during or after the electroweak 
transition. 

Without the bound (|TT|) on the singlet majoron model, it would be natural to assume 
that v s and M are large. For T BL < M ~ v s , the lepton- violating decays and inverse decays 
would have gone out of equilibrium before baryogenesis, evading the bound ([14]). Then one 
would be left with only the weaker constraint coming from the dimension-5 operator fllTf). 
Until a nonvanishing neutrino mass is measured, this puts no restriction on the scale of 
baryogenesis. But with the restriction (O) the story changes dramatically, and the window 
of opportunity for producing the B — L asymmetry becomes very narrow. With (|Tl"| ) the 
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bound from decays and inverse decays is the only one available as the dimension-5 operators 
go out of equilibrium at T ~ 10 TeV, for m u ~ 25 eV, and cannot be used at T > v s . 

The original derivation of the bound (fLip was incomplete in the sense that the authors 
did not solve, nor even formulate the Boltzmann equations for the asymmetries. Because 
we are studying a theory with spontaneous symmetry breaking, where the lepton-violating 
decays leading to ( |14D only exist below the critical temperature T c at which v s starts to 
develop, and since we are finding that v s must be near the weak scale, it behooves us to 
reexamine in greater detail under what conditions the validity of flli]) is generally assured. 

Our basic assumption is that all the other processes proceed fast compared to those 
that violate lepton number, so that at each instant of time an equilibrium between various 
chemical potentials is attained. The evolution of all asymmetries is then controlled by the 
evolution of the individual lepton asymmetries due to the decay processes N <-> HL. 

It is convenient to choose the mass eigenstates of the right-handed sector as the two 
helicity states of the previously defined Majorana field iV = vr+vr c ., which does not mix with 
the light states until temperatures below the weak scale. Above the critical temperature of 
the phase transition in the singlet sector these are precisely the eigenstates of the conserved 
right-handed neutrino lepton number. Moreover it is convenient to work in the basis in 
which the Majorana mass matrix is diagonal. Since we are interested in temperatures above 
the weak breaking scale, all the other fermions remain massless. As the weak symmetry 
eventually gets broken, the light neutrinos gain a seesaw mass matrix whose diagonalized 
form is 

m v = U T m T —mU, (15) 
where M is the heavy Majorana mass matrix, m is the Dirac mass matrix 

m = hv/y/2, (16) 

and U is the neutrino mixing matrix that relates the flavor and mass eigenstates. 

The Boltzmann equation for the evolution of the vg number density, due to the lepton 
number violating decays is given by 



S 



(d t -3H)n Ut = J Y[j^ F (27rmp 1 -p 2 -p 3 )^:\M\ 2 1 

i=l \ J i j 

x fh(p2,^h)fu(p3,^)e m 

x { /i+ (^,^.)( e -^- e ^-^)) (17) 

where £ is the generation index. The first term in the curly brackets in ( |T8"| ) comes from 
iVj + and the second from Nj- decays, with the index j running over all heavy neutrino 
species and H is the Hubble expansion rate of the Universe. The maximal lepton number 
nonconservation below T c in the decays N± <-> Hvl-, H*vl c is reflected in the fact that the 
helicity states N± decay with equal total rate to both allowed final states with opposite 
lepton number. Indeed, the matrix element for the decay of Nj\ is given by 

\M\ 2 jX = {hjtfMjfa ■ V v ± MjSj\ ■ Pv ), (18) 

where Sj\ is the helicity four-vector of Nj and — (+) refers to the Hpl c (H*ul) final state. 
The various decay channels differ by the spin dependent part in (|18|) , but it is easy to see 
that this contribution identically vanishes upon the phase space integration. 

Subtracting the Boltzmann equation for the number density of vfs from that of z//s and 
linearizing in the chemical potentials, one finds 

(d t - ?>H){n n - n ve ) = -4 — - — x ^ h H h u I u . (19) 

The function In can be evaluated exactly to yield 



HI 



oc 



COsh(tt£M + ~j\Ju 2 — 1) 

1 + e ux ) 2 \cosh.(ahU — / yy/u 2 — V 



du— r- In 



x smh(a e u + jVu 2 - 1) j ^ 
sinh(a^w — 7v / w 2 — 1) j 



where a e = {M 2 +m 2 Ue ~m 2 h ) / '4M;T ', a h = (Mf +m 2 h — m 2 Vi ) / 4M;T and 7 = \ 1/2 (M 2 , m\, m 2 h ) 
/4MjT, with A(:r, y, z) = (x — y — z) 2 — Ayz. This is a generalization of the inverse-decay 
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rate used in J29|. Note that the Ni chemical potentials do not appear in the equation fll9p. 
This is again due to maximal lepton number violation in the decays; should there be a 
difference in the decay rates, the Ni asymmetries would be nontrivially coupled to zVs. In 
the present case the evolution of the Ni asymmetry decouples and the Ni asymmetries are 
driven to zero independently of (and faster than) the asymmetries of the v^s. A subtlety in 
the finite temperature calculation of the decay rate that deserves to be mentioned is that the 
temperature dependent masses enter only through the kinematics and the dispersion relation. 
In particular, the spinor functions must be taken identical to the vacuum spinors, but with 
a modified dispersion relation [p8| . In all our calculations we have used the approximate 



dispersion relation E = yp 2 + m 2 (T) where m(T) is the appropriate temperature dependent 
mass. 

Equation (|20|) can be greatly simplified if one substitutes the Maxwell-Boltzmann ap- 
proximation for the phase space distributions. In general this is not legitimate, because 
then, for small Mi/T, the enhancement factor due to the singularity in the boson distribu- 
tion function at E^/T ~ would be ignored In practice however, this singularity is shielded 
by the temperature dependent mass of the final state particles and the MB approximation 
becomes effectively very accurate. We then find the approximate expression 

hi Wf *i [jr) , (21) 

where K\(x) is the modified Bessel function. This result coincides with (|20| ) to within few 
percent over the whole range of interest, and becomes exact for large Mi/T. 

The Higgs chemical potential fih can be solved in terms of the leptonic chemical potentials 
from the equilibrium conditions imposed by the fast weak processes. Following P], and using 
fJ>t = Hu e due to fast weak interactions, we find 

MA=^-f>- ( 22 ) 

Inserting this result into the equation (|l9|) and expressing the chemical potentials in terms 
of the asymmetries in the comoving volume, we obtain 

K = n -^^ ^ (23) 
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where s is the entropy density and g* is the number of relativistic degrees of freedom. We 
find the following coupled set of equations for the lepton asymmetries 

with 

24 2 (M^+m|-7^ 1/2 2 2 2 /MA 

Eq. (p4|) is valid as such for both neutrino and left handed charged lepton asymmetries, 



which is why we replaced the subscript vi with £ in ( f24|) and fl25|) . One should note that (|24|) 
allows for a nontrivial equilibrium solution (some of the L/s nonzero) only if at least one 
of the rates Ti vanishes, or is small enough fl3| . In the opposite case all L/s and hence all 
the other asymmetries as well are driven to zero with a rate controlled by (f2~4"|). The baryon 
asymmetry is related to the sum of the individual neutrino or left handed charged lepton 
asymmetries by 

B{t) = - A -Y, L ^)- (26) 
6 t 

Given the initial asymmetries and the various couplings that appear in the expressions for T^, 
one could compute numerically the baryon asymmetry as a function of time from equation 
(|24|). For the purpose of obtaining a bound on the neutrino mass, we can make the following, 
most conservative, simplifying assumption about the parameters involved in ([24]): "the initial 
asymmetry corresponding to the smallest of rates is the largest." One should bear in mind 
that the fast electroweak processes enforce the initial conditions for the slow processes such 
that all initial asymmetries are of the order as L/s. The faster decaying asymmetries in 
eq. fl24|) rapidly reach their fixed point values and it can be integrated exactly for the slowest 
one to yield: 

/ 33 /"* \ 
L e (t) = L e {t ) exp (-— J dt T t (t)J . (27) 

The bound on the neutrino mass arises from the requirement that the present baryon asym- 
metry exceeds the known lower limit of B m i n ~ 4 x 10~ n , which in terms of the lepton 
asymmetries becomes 



3 



4 

> ~B mm , (21 
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where the time tEw refers to the time of electroweak phase transition. With our conservative 
assumption we then obtain the bound 

"" AtTl ) slbm^ir, (29) 



t 



mill 



33 V 4-10-11 , 

using in the last step a reasonable maximum value of L° imax ~ en 7 /s £j 10~ 2 for the initial 
asymmetry, where e is some generic CP-asymmetry parameter f2~9|| . Up to now our discussion 
has been quite general and eq. (|29|) is valid for either explicit or spontaneous symmetry 
breaking. Assuming that the temperature dependent mass of the N{ is given by 

M 2 (T)=M 2 0l (l-(T/T c ) 2 ) } (30) 

we can write the integral appearing in ( |29| ) in the spontaneous symmetry breaking 

J dt T e ~ 284 m* tt —mu x I (31) 



where 

d^ T — — /W^iW, (32) 

with z = Mi(T)/T, f(z) = (1 + rf - r^AV^l, rf , r 2 )#(l - r £ - r A ) and r a = m a (T)/M t (T). 
The explicit symmetry breaking case can be obtained by taking the limit t c — > 0, T c — > oo. 
In the singlet majoron model the various temperature dependent masses that we use in our 
calculations are given by 

M 2 (T) = 2h 2 2l v 2 s {T) ; v 2 (T) = v 2 (l - T 2 /T 2 ) 

mj(T) = ±(3g 2 + g' 2 )T 2 (33) 

m 2 h {T) = l(3g 2 + g' 2 + Ah 2 + 8X H + ^X sH )T 2 -X H v 2 -^X sH (v 2 -v 2 (T)), 

where g and g' are the usual weak and U(l)y couplings, h t is the top quark Yukawa coupling, 
Xh, X s , X s h, v and v s are the quartic scalar couplings and VEV's defined in the equation (||) 
and T c is the (B — L)-symmetry breaking temperature given by 

V 2X s + X SH + h 2 2 J 
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where h 2 = J2j The expression for m/ l (T) in (|34| ) assumes that for temperatures T c > 
T > Tew, the field is located in the temperature dependent asymmetric miminum in the 
singlet direction with the electroweak symmetry still unbroken, (\a\) = v s (T)/2 and (\H\) = 
0. We have computed the integral / for a wide range of parameters and find that it is 
essentially constant I ~ 4, for Mj > m^. The constant behaviour arises because the integrand 
~ z 3 Ki(z) peaks for rather large values of z ~ 3 so that the effects of final state masses and 
the Moi/T c dependence are strongly suppressed. Thus the final bound arising from fl29f ) is 
essentially the same for either spontaneous or explicit breaking. With a constant / the LHS 
of the equation (|29|) becomes proportional to the quantity [rwM~ x, m)a- 



Our main result now comes from combining ( p8|) and fl29|) , using the constant value of I 
just determined: 

nun {m ] M- l m)u < 8 x l(r 4 ln ( - ^°_ n ) e V ^ 1<T 2 eV , (35) 

where again, the maximum value for L was used in the last step. Notice that the careful 
treatment leading to eq. ([35]) has yielded a somewhat less restrictive bound on the scale of 
the seesaw neutrino masses, m 2 /M, than previous investigations. But another interesting 
feature of (|35|) is that, unlike previous derivations, it does not give any direct bound on the 
masses of the neutrinos, because m)M~ l m is different from the mass matrix m T M~ l m. For 
example, consider just two generations and assume that M~ 1//2 m is proportional to the Pauli 
matrices <j\ + a 2 - Then m^M _1 m is diagonal and nonzero, but the mass matrix vanishes 
identically, and thus the bound (|35f) is not directly related to the masses. Such a situation 
could arise either due to fine-tuning of the mass matrices or (approximate) lepton flavor 
symmetries. Otherwise one generically expects m) M~ l m to be related to the masses and 
mixings via 

3 

(mtM- 1 ™)** ~ X>jEk| 2 , (36) 



as can be seen by writing M 1//2 m = V^DU, where D is diagonal, and assuming that V 
is either close to the identity, or has off-diagonal elements of the same order as those of U. 
Then m)M~ l m = W\D\ 2 U, whereas the mass matrix is U T DV T VDU. If V T V = 1 + e, 
then (pop is correct to lowest order in e. If e is large but the off-diagonal elements of U are 



13 



also large, then (^) will again be approximately true, since it is dominated by the largest 
eigenvalue of D, that is, the heaviest neutrino mass, regardless of the generation index I. 

If either of the reasonable assumptions leading to eq. fl36|) holds, we can combine the 
baryon-preservation bound (|35| ) with the majoron decay bound (|Tl|) to obtain what is po- 
tentially a much more stringent limit on the majoron scale than was given previously in 
eq. (|il~l). Eq. (|35|) tells us that the mass of the heaviest neutrino Vh into which the ma- 
joron can decay should satisfy m Uh < 10~ 2 U^ eV, where £ is the neutrino flavor whose 
conservation is most weakly violated. Then (|TT|) becomes 



v < 



max (200 GeV(27Y)~ 4/7 £/ w 8/ V 1/7 , vj . (37) 



Note that we cannot obtain a bound for v s < v, since sphaleron interactions will have 
decoupled by then. 

A further application of our bound is a constraint on the Yukawa couplings of the right- 
handed neutrinos to the Higgs doublet, about which no information has hitherto been avail- 
able. Rewriting the mass matrices in terms of Yukawa coupling matrices and using the 



constraint (|37D on v s we obtain 

J2 < 6 ■ 10- 14 max (o.8(27Y)- i/7 U he 8/7 X' 1 ' 7 , l) . (3* 



mm 



Even if the gravitationally induced effects (parametrized by A) or the mixing Um are very 
small, this is quite a strong result, for it says that either all the Yukawa couplings ha are 
even smaller than that of the electron or there exist large hierarchies within hu with some 
of the elements being extremely small; hu 10 . 

The dramatic results ( |37| , |3~8|) apply in the case of generic mass matrices for the neutrinos, 
but in model building it is often interesting to consider the nongeneric case where some 
combination of flavors L is approximately conserved, for example L = L e — + L T . If L 
is only weakly violated, then the decay vertices we have considered have only very slow L- 
violating channels, since the heavy decaying neutrinos iV have pseudo-Dirac masses of order 
M (which almost conserve L) rather than Majorana masses that would violate L maximally. 
The decays are therefore predominantly of the form iV — ► Hu and iV — > H*v. The fastest L 
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violating process are N-N oscillations coming from the small Majorana mass terms fi |3Q| 



followed by the L-conserving decays. These are faster than any oscillations ocurring in the 
light neutrino sector because the mass splittings, hence the oscillation rate, are much larger 
for the heavy neutrinos. To be exact, we would have to write the Boltzmann equation also 
for the asymmetry in N and combine it with that for v in order to track the asymmetry in 
L. But roughly, we can account for this effect by multiplying the previously computed decay 
probability by the probability of an oscillation, given by 

P = sin 2 (2^sin 2 (|§) (39) 

where for pseudo-Dirac neutrinos the mixing is maximal (sin 2 (2$) = 1), and Y is the rate of 
interactions of N which will damp the oscillations if they are much faster than the oscillation 
rate SM 2 /2E = \xMjE. T is given by the sum of rates J2e ^ which we computed previously, 
eq. For very fast oscillations the second sin 2 () factor averages to 1/2, but if \i is small 

then our bound (B3) starts to disappear. The condition for Q3"5p to be valid in the case of an 



approximate symmetry is therefore that \iM ,> TT at temperatures of order M, or in terms 
of the dimensionless paramter \xjM that quantifies the small departures from L symmmetry, 

^ > ^EIM 2 -^^. (40) 

Therefore even in the extreme case of M — 100 GeV and m v . A = 10 MeV, a very small 
symmetry breaking fi/M ~ 10~ 6 would be sufficient for our results to still apply. 

5 Caveats 



Finally, we discuss some special situations, where our bounds 037|) and fl38|) could be avoided. 
One such situation is the case in which the baryon asymmetry is actually generated below 
T c . Since the wash-out due to the decays is not effective for T <C Mj(T), the asymmetry 
generated at Tbl is safe if Tbl *C T c . However, due to the bound (|TT|) on v s with T c ~ v s 
this should be seen as an extremely tight constraint on the scale Tbi)- 

Another possibility is that the baryon asymmetry is generated at some high scale, and 
the electroweak symmetry breaks before the B — L symmetry even though v s > v. If this 
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was the case, then there never was a period where the sphaleron interactions and the lepton 
number violation were in equilibrium together and the bound fl36|) would not be valid. To see 
how generic this is, we must compare the critical temperature for the electroweak breaking 
(at the field origin) with the (B — L)-symmetry breaking scale T c (0). We find 




where r = (2M^ + M| + 2mf)/v 2 . Note that if the B — L symmetry is indeed broken first, 



from the origin after (B — L)-breaking. The parameters of the quartic couplings in the scalar 
potential must always satisfy the relation X s h < 4AA S in order for both symmetries be broken; 
otherwise there would be no light neutrino masses at low scales. Moreover, the requirement 
that the one-loop corrections do not render our vacuum unstable introduces a bound on the 
Yukawa couplings h 2 i. This bound can be naively estimated by computing the coefficient of 
the \a\ 4 In \a\ term and requiring that this coefficient be nonnegative: (§Ag+ ^A s #— J2j ^fj) > 
0. One can then show that the aforementioned constraints with v < v s and Tew > T c would 
be satisfied only in a small region of parameter space, where A s <C X s h <A~0(1). 

There are two more subtle cases where the bound fl35|) could be weakened, related to the 
computation of the integral I in (|32|). First, it was assumed that the parameter M Qi /T c is not 
excessively large. Using the vacuum stability bound and the bound on the quartic couplings, 
one can show that M 0i /T c >, some number a would require y/3a(l + (x") 1//2 )~ 1 ^ 2 ~ h 2 ^ 2X S . 
This inequality has solutions for a ^ 1, i.e., only for nonperturbative values of couplings 
X,\ s ,ti2 ^ 0(1). A more severe problem arises if Moi/T c is excessively small. As noted 
earlier, the integral in (|32p gets its largest contribution from the region z = Moi(T)/T ~ 3. 
For small M 0i /T c this peaking occurs at temperature T p ~ M 0i /3. It is always possible 
to make M 0i much less than T c merely by letting h 2 become small. When M 0i < 
our bound scales roughly as h^ 1 . However, when v s is small, a very small value of h 2 is 
rather unattractive anyway, since the mass becomes very small, rendering the seesaw 
mechanism less natural. 





(41) 



the temperature ([TI]) is not the true electroweak breaking temperature, since the field shifts 
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6 Conclusions 



While it is clear that small neutrino masses are an interesting and potentially necessary 
extension to the standard model, the origin of these masses remains unclear. Many of the 
possibilities involve lepton number violating Majorana masses, as is the case in the singlet 
majoron model discussed here. The combination of lepton number violation in conjunction 
with the B + L violation in the standard model leads to the grave possibility of washing 
out any primordial baryon asymmetry. Thus, the survival of the baryon asymmetry enables 
one to constrain models such as the singlet majoron model (unless, of course, the baryon 
asymmetry is produced very late). We were thus able to strengthen the bounds of ref. |17j . 
We have shown under general circumstances that at least one neutrino mass is constrained 
to have a value m v < 10~ 2 eV and the neutrino Yukawa coupling is constrained to h < few x 
10~ 7 . The singlet VEV is correspondingly constrained by v s £ v. While the majoron model 
is still viable with a low value of v s , models with a large VEV must be carefully tuned so as 



to avoid the constraints discussed here as well as in IT7 
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